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Interaction is so ubiquitous that imaging a 
world free from it is a difficult fantasy exercise. 
At the same time, in understanding any com¬ 
plex physical system, our ability of accounting 
for the mutual interaction of its constituents is 
often insufficient when not the restraining fac¬ 
tor. Many strategies have been devised to con¬ 
trol particle-particle interaction and explore the 
diverse regimes, from weak to strong interac¬ 
tion. Beautiful examples of these achievements 
are the experiments on Bose condensates mm , or 
the recent experiments on the dynamics of spin 
chains [4][5]- Here I introduce another possibility, 
namely replacing the particle-particle interaction 
with an external stochastic field, and once again 
reducing the dynamics of a many-body system to 
the dynamics of single-particle systems. The the¬ 
ory is exact, in the sense that no approximations 
are introduced in decoupling the many-body sys¬ 
tem in its non-interacting sub-parts. Moreover, 
the equations of motion are linear, and no un¬ 
known external potential is inserted. 

The idea of replacing the many-body system under in¬ 
vestigation with a non-interacting doppelganger is not 
new. From a theoretical point of view, a starting idea 
has been to treat the interaction as an external pertur¬ 
bation. Interaction “dresses” the particles and new funda¬ 
mental particles appear for the description of a physical 
phenomenon. This is the tenet of the Landau’s theory 
of the Fermi gas mapping a strongly interacting electron 
gas into a system of weakly interacting quasi-particles 
[3 0. More modern approaches replace the particle- 
particle interaction with an external effective potential, 
e.g., the Thomas-Fermi’s theory and the Hartree-Fock 
approximation jS] which have all converged now some¬ 
what in the Density Functional Theory 13 IH HH1- The 
price to pay for this huge simplification is the inclusion 
of a unknown non-linear potential in the dynamics of the 
fictitious non-interacting system BE]. Density func¬ 
tional theory has been instrumental in understanding 
many physical, chemical, and biological phenomena at 
the nano-scale and in augmenting the theoretical predic¬ 
tion potential. 

With hindsight the results I will present in the fol¬ 
lowing are not completely surprising: for example when 
dealing with magnetic systems, a common approximation 
consists in replacing the dynamics of spin operators with 


the dynamics of their quantum averages. Often, these 
averages have a random behavior since they mostly con¬ 
sist of the superposition of a static magnetic moment and 
small dynamical fluctuations. 

In this Letter, I will show that this analogy is even 
more stringent and can be made exact for the case in 
which particles interact via a potential that depends on 
one operator of particle i multiplied by an operator of 
particle j. Indeed, the dynamics of such a system is 
exactly equivalent to the dynamics of a system of non¬ 
interacting particles in the presence of known stochastic 
potentials. The dynamics of the many-body system is 
then recovered by constructing the proper wave-function 
or density matrix and averaging the results over the 
stochastic fields. This study finds direct application in 
the dynamics of arbitrary spin-chains and gives immedi¬ 
ate access to the exact high-order correlation functions 
that recently have been probed experimentally m 0- 
Moreover, the results presented here open up the pos¬ 
sibility of investigating the dynamics of large systems. 
It is well known indeed that the numerical solution of 
the equation of motion of the many-body density ma¬ 
trix, p, is limited by its intrinsic large dimensionality 
which scales at least quadratically with the number of 
states one needs to consider, compared with the linear 
scaling of the wave-function. Our result shows that for 
certain cases, one can bring the exact dynamics down 
to the evaluation of N relatively small density matrices, 
thus regaining a linear scaling. To give an estimate of the 
problem with the interacting system, with a chain of N 
spin 1/2, the exact many-body states belongs to a space 
of dimensionality 2 N , and therefore the density matrix 
has dimension 2 N x2 N . It should be then apparent that 
we can investigate numerically small chains, routinely up 
to N ~ 10 — 15, after which the computer memory re¬ 
quirements will be prevailing. Our result brings down 
this requirement to store up to N 2 x 2 matrices. 

Let us begin with considering a many-body quantum 
system, whose dynamics is determined by the many-body 
Hamiltonian 

N N 

H = ^2 hi + X ® Xj, (1) 

i i<j 

where hi is a single particle Hamiltonian and Xi is some 
operator acting on the particle i, Aij the interaction con¬ 
stant between particle i and particle j, and N the total 
particle number. By Newton’s third law, the interaction 
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constant is symmetric, i.e, A ij = \j t i. The many-body 
density matrix of the system evolves according to the von 
Neumann’s equation 


hd t p(t) 


— l 




( 2 ) 


For simplicity, let us assume that the initial density ma¬ 
trix is the direct product of single particle density matri¬ 
ces 


Pit = 0) = 0 Pi(t = 0), (3) 

i=l,N 


then we can find a set of N 2 independent white complex 
noises cotj ( t ) for which the single particle density matrix 
Pi evolves according to m 


Mpi(t ) = - i [hi, f>i(t)\ dt 

+ ^2 ^ 9 ([»», Pi] duiij - i {cti, pi} dui*i 

3=1 


(4) 


and the exact total density matrix is given by 

P0) = 0 Pi(t) (5) 

i=l,JV 

where the 7TT denotes the average over all the white 
noises. In Eqs. ft and ([ 4 ]), [A, B] = AB — BA and 
{A, B} = AB + BA are the standard commutator and 
anti-commutator of any two operators, respectively. The 
white noises u>ij are complex Wiener processes that sat¬ 
isfy 


dUijdUk,p — 2 di^Sjpdt, ( 6 ) 

where = 1 if i = k or vanishes otherwise mm- 
We prove Eq. (J5| starting from Eqs. Q in the Methods 
section. 

Let me now discuss two important points. First, it 
may appear that the initial correlation between the par¬ 
ticles is lost in this theory. This is not the case. In¬ 
deed, the following discussion can be generalized to the 
case in which the initial condition is given by p(0) = 
® i=i n 7 &Pi (0)> where 7 k is a set of coefficients such 
that = 1, to ensure that if tr/3*(0) = 1 for any 

i, then trp(0) = 1. The total density matrix in this 
case, owing to the linearity of the equation of motion, 
will then be given by p(t) = Efc 7fc0 i= i,jv Pi(t)i where 
each p\ evolves according to Eq. Q with initial condition 
P?(0). For simplicity, in the following I will maintain that 
Eq. (§ is satisfied by the initial density matrices. The 
second point is the form of the particle-particle inter¬ 
action: if any operator of particle i commutes, or anti¬ 
commutes, with any operator of particle j then one can 
always expand any two-particle interaction as a series of 


products of single particle operators like in Eq. ([T]) . The 
theory we are putting forward then is useful for those 
case in which the particle-particle interaction is written 
as a finite sum of products of pairs of single-particle op¬ 
erators. Finally, it should be clear that even if the initial 
density matrix pi(0) describes a real particle, its time evo¬ 
lution pi(t) cannot be associated with the dynamics of a 
real particle. This is easily seen by the fact that Eq. Q 
does not preserve either the positivity or the unitarity of 
Pi{t), even if we assume pi(0) is a definite positive matrix 
of unitary trace. We will discuss later on how to reduce 
the many-body density matrix to a proper single-particle 
(reduced) density matrix that can be used to investigate 
the single-particle properties of the system. 

The exact dynamics of the many-body density matrix 
usually contains a redundant amount of information. For 
practical purposes, it is usually more convenient to trace 
out some of the degrees of freedom and obtain the expec¬ 
tation value of single- or two-particle operators. Within 
this theory, this procedure emerges naturally from the 
definition of the operator pi. We have for example for 
the time evolution of the subsystem j, 


Pf{t) = Pj(t)Y[tTpi(t), (7) 

i-Aj 

where tr indicates the trace operation on the matrix pi. 
Notice that in general tr pi(t) is a function of time, and 
we cannot expect that it equals 1 at all times. On the 
other hand, we expect that tr pf{t) = tr p(t) = 1 at all 
times, so that pj* does define a proper density matrix for 
the subsystem j m 

As a first example of application we consider the case 
of two interacting spins. We assume there is not any 
external magnetic field. The Hamiltonian for this simple 
system is H = Xaf g> erf where of is the 2x2 third Pauli 
matrix. For simplicity and without any loss of generality 
we can set A = h. According to Eq. Q we need to solve 
the two independent stochastic master equations 

dpi = - [a z ■ pi(dui2 — idc*4i) — pi ■ a z {idu>2i + duji 2 )}, 

dp2 = 2 \ aZ ' P2{dt02i — idhj\ 2 ) ~ P‘2 ' cr z {idu>\ 2 + duj 2 1 )] • 

( 8 ) 


One can easily prove that the solution of these two 
stochastic equations are 


Pi(t) 


e -u>12p21 


e ui 12 pl2 \ 

e ~ iu] 2i p22 J ■ 


(9) 


Here, we have introduced explicitly the elements of the 
initial state p‘\ l . / 02 (f) can be obtained from this ex¬ 

pression by substituting 1 —> 2 in the subscript of 
the initial state, and by swapping 1 2 in the sub¬ 

script of to. Some straightforward algebra now leads 
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to the expression of p(t), and with the properties of 
the averages we will discuss in a moment, to the fi¬ 
nal expression for the total density matrix p(t). We 
can prove that the density matrix obtained in this way 
is identical to the one obtained by evaluating p[t) = 
exp(— iHt)p(Q) exp(iUt) where p(0) = p\(0) 0 /S 2 (0). To 
obtain these results, we need to evaluate the stochastic 
averages of two functions: namely, both exp(wi 2 +W 21 ) 
and exp(wi 2 + wu* 2 ). To do that, we observe that for any 
real Wiener process, p , and any complex number, a, we 
have exp(ctp) = exp(— a 2 t/2) [T5] . We easily then arrive 
at exp(wi 2 + W 21 ) = 1 and exp(wi 2 + zw* 2 ) = exp(—2 it). 
We can now evaluate the reduced density matrix for the 
spin 1, pf (t). Again with some straightforward algebra, 
and assuming that trpj/O) = 1, we obtain 

Pl() V(p1 2 w)* pf J (0) 

where p\ 2 (t) = p\ 2 (cos 2t + i(pf — pf) sin 2t). With 
the reduced density matrix we can therefore calculate 
the expectation value of any single spin observable, 
e-g (erf) = pf - pf and (erf) = 2Re(pf 2 ) cos(2t) - 
2Im(p{ 2 )(p 22 — pf 1 )sin(2t). 

The generalization to a chain of N interacting spins 1/2 
therefore follows in the same footsteps. If the start with 
the Hamiltonian H = 1 Aj 0 "/ 0 erf, with Xjj = 

0 and use Eq. Q, we obtain the single spin stochastic 
density matrix 


Pi(t) = 



( 11 ) 

From this result, deriving the single spin reduced density 
matrix is now straightforward, 


Pi (t) = n 1 *Pk{t)pi{t) = ^ j P ^ t y ^ , (12) 

k^i 

where we have introduced the functions f t (t) = 

Pi 2 nf=i, k^i ( Pk 1 exp(-2«A ik t) + pf exp(2i\ ik t)). In 
the same way, we can build the 2-spin reduced den¬ 
sity matrix starting from its definition pf/ (t) = 

n^Li tr p k (t)pi(t) <3 Pj(t). The knowledge of pfft) al- 
k^i,j 

lows us to calculate the correlation functions that have 
been recently measured mm, in accordance to new the¬ 
oretical results mm- The present theory provides an 
easy way to reproduce and generalize those results. 

We can investigate how the presence of a magnetic 
field, for example in the x direction affects the dynamics. 
To do that, we consider two interacting spins with the 
total Hamiltonian given by H = h\af + /i2<7 2 + Acrf 0 erf. 
We can write down the stochastic equations of motion 
for the two single-spin density matrices, but their ana¬ 
lytic solution of little use since a z and a x do not com¬ 
mute and we have to revert to a numerical solution of the 


equations of motion. In Fig. [T] we report the dynamics 
of a few elements of the density matrix pf (t) calculated 
using Eqs. |7| and Q, after taking the average of 10 6 
independent realizations of the white noises. As initial 
condition we have assumed the two spins are in the mixed 
state, pf = pf = 1/2 for any i and j. 

The possibility to measure in time, and locally the time 
evolution of the correlation between two spins has re¬ 
cently captured a lot of attention. HG3 In the experi¬ 
mental set-up, we could change the spin-spin interaction 
in such a way to explore the transition between a XY- 
model to Ising model. In particular, this possibility has 
been investigated in linear spin chains. Of particular in¬ 
terest is the spin-spin correlation function, since it can 
be seen as a way to measure the spin propagation speed 
and correlation time. To show how the present formalism 
can be applied to this case, we consider the Hamiltonian 
H = hi<r? + 1/2 ( a i 0 describing a 

chain of N interacting spins in the presence of a static 
magnetic field Bi = (0,0, hf. In the following, we con¬ 
sider the case of a chain of 11 spins, in the presence of a 
uniform magnetic field, hi = h for any i, and where the 
interaction if restricted to first neighbors. A quantity of 
interest for these systems is the time evolution of the cor¬ 
relation function, defined as Cij = (plof) — (of)(crJ). In 
Fig. we report the time evolution of Cej(t). We have 
chosen as initial condition that all the spin but 6 are in a 
mixed state, while at t = 0 the state of spin 6 is up, that 
is p e ’ = 1, p 6 ' = p 6 ’ — p 6 ’ = 0, while we have used 
h = 1 and A,y = 0.01 if |i — j\ = 1, or 0 otherwise. We 
have chosen a time step of At = 0.002 and averaged over 
10000 realizations of the stochastic noise. It is seen that 
the correlation grows rapidly for the spin 5 and 7, while 
for the other spins, that are not directly connected with 
spin 6, it remains rather small. This can be compared 
with the results of the experiments reported in 110 . 

It is well known that the von Neumann equation is 
equivalent to the Schrodinger equation to almost any pur¬ 
pose in standard quantum mechanics 0 H3. It then 
appears natural to extend the result of this Letter to 
the many-body problem described by the time evolution 
of the many-body wave-function, |’F(f)). We can prove, 
following essentially in the footstep of the proof given 
for the many-body density matrix, than it is possible 
to find stochastic single-particle wave-functions \ipi) such 
that, if at the initial time |v&(0)) = (S)iL~i l'0i(O))i then for 

any subsequent time |’F(t)) = IV’i(^))- The states 

fi(t) evolve according to a stochastic Schrodinger equa¬ 
tion mimg which has a form similar to Eq. Q. Inter¬ 
estingly, while the standard relation p(t) = |’F(t))(\E'(t)| 
is satisfied, there not exists a similar relation between 
the single-particles, i.e., in general we should expect 
that pi(t) f \ijji(t))(il>i(t)\. Indeed, one can prove that 
\ipi{t ))has an equation of motion that does not re¬ 
duce to Eq. Q. Finally, we would like to comment that 
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a formalism based on the wave-function | tpi) seems lim¬ 
ited, since it is not clear how one could possibly obtain 
information on the single-particle properties of the real 
many-body problem. In fact, if we start with the defi¬ 
nition of pf (t) from Eq. Q, and use the wave-function 
|4/(t)) instead we get, 

pf(t) = trp(i)| but j = tr(|tf(i))(tf(i)|) | but 3 

( ~N N 

0 hM*)) 0(^0)! 

i= 1 i= 1 

It is therefore not possible to swap the trace operation 
with the average over the realizations of the stochastic 
noise, therefore precluding an alternative way to Eq. (J7|. 
Similar problems arise when one considers any ro-body re¬ 
duced density matrix as obtained from the single-particle 
stochastic wave-functions. 

In conclusion, we have shown that the dynamics 
of a many-body system, with multiplicative two-body 
particle-particle interaction can be reduced to the inves¬ 
tigation of the dynamics of N single particle stochas¬ 
tic systems. We have shown how to calculate any re¬ 
duced ro-body properties starting from the solution of 
these stochastic dynamical equations. We have applied 
this formalism to the case of spin chains in the presence 
of a finite magnetic field. 



METHODS 


Here we present the proof of Eq. ([5]) , based on the Ito 
calculus mum. 

Proof. The proof of Eq. ([5]), and that p(t) follows 
the dynamics induced by the von Neumann’s equation 
([2]) is a straightforward generalization of the results of 
Shao on the dynamics of a quantum system in contact 
with an external environment EH HU. Let us consider 
the stochastic total density matrix, p(t) = (^) i Pi{t). Ac¬ 
cording to Ito’s stochastic calculus, the dynamics of p is 
determined by, dp(t) = J2jLi Pip ® dpj <g> p 3+ pjv+i + 

J2k<j=i (Pi,k ® dp k <g> pk+ij ® dpj (g) p j+ i,jv+i) , where 
we have introduced the short-hand notation p i3 = 
The proof that p(t) = p(t) then continues by 
using Eq. Q, Eq. |6]), and [x l) p 3 \ = 0 if i ^ j, to arrive 


at dp(t) 


H,p(t) 


dt. Due to the linearity of this 


equation of motion, and that by definition p(0) = p(0), 
we conclude that p(t) = p(t) at any time t > 0.D 
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FIG. 1. The time dynamics of the real parts of (pf)n and 
(pf) 21 calculated from Eqs. <|t]) and ([ 4 ]). As initial condition 
we have set the two spin in the “mixed” state. We have av¬ 
eraged over 10 6 independent realizations of the white noises. 
For the calculation we have chosen hi = h 2 = 1 and A = 0.02, 
and a time step, At = 0.01. We can see that after t ~ 100 the 
finite number of realizations and the finite numerical accuracy 
in solving the stochastic differential equations are introducing 
an error in both p} 1 and p\ 2 . 
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FIG. 2. Dynamics of the correlation function Ce,j(t) in time. 
We have chosen for the parameters h = 1 and A = 0.01. All 
the spin but 6 are in a mixed state, while spin 6 starts in the 
state “up”. We have chosen a time step of At = 0.002 and 
averaged over 10000 realizations of the stochastic noise. 






